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Abstract 

Purpose: Many planning methods for high dose rate (HDR) brachytherapy treatment planning require 
an iterative approach. A set of computational parameters are hypothesized that will give a dose plan that 
meets dosimetric criteria. A dose plan is computed using these parameters, and if any dosimetric criteria 
are not met, the process is iterated until a suitable dose plan is found. In this way, the dose distribution is 
controlled by abstract parameters. The purpose of this study is to improve HDR brachytherapy planning 
by developing a new approach that directly optimizes the dose distribution based on dosimetric criteria. 
Method: We develop Inverse Planning by Integer Program (IPIP), an optimization model for computing 
HDR brachytherapy dose plans and a fast heuristic for it. We used our heuristic to compute dose plans 
for 20 anonymized prostate cancer patient image data sets from our clinic database. Dosimetry was 
evaluated and compared to dosimetric criteria. 

Results: Dose plans computed from IPIP satisfied all given dosimetric criteria for the target and healthy 
tissue after a single iteration. The average target coverage was 95%. The average computation time for 
IPIP was 30.1 seconds on a Intcl(R) Core™2 Duo CPU 1.67 GHz processor with 3 Gib RAM. 
Conclusion: IPIP is an HDR brachytherapy planning system that directly incorporates dosimetric crite- 
ria. We have demonstrated that IPIP has clinically acceptable performance for the prostate cases and 
dosimetric criteria used in this study, both in terms of dosimetry and runtime. Further study is required 
to determine if IPIP performs well for a more general group of patients and dosimetric criteria, including 
other cancer sites such as GYN. 



1 Introduction 

High dose rate (HDR) brachytherapy is a radiation therapy for cancer. In HDR brachytherapy, radiation is 
delivered directly to the tumor site via a moving radioactive source in temporarily inserted catheters. Dose 
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4 is controlled by altering the dwell times, the time spent at points along the catheters. Studies have shown 

5 that brachytherapy is a highly effective treatment [TJ[21[31IU|5Jini[71[E]- 

e Inverse planning explicitly uses anatomical information and dwell positions when computing HDR brachyther- 

7 apy dose plans. Clinical experience with HDR brachytherapy has led to organ doses that are correlated with 

s biologically acceptable results [9! , more specifically, the eradication of the tumor with reasonable side effects 

§ to healthy tissue. Recently, these organ doses have become the dosimetric criteria that form the objective of 

10 many dose planning systems. 

11 Many planning systems provide tools to quickly evaluate dose distributions and guide the user toward 

12 a final dose plan. However, finding a suitable dose plan may take several attempts. During dose planning, 

13 a set of computational parameters are hypothesized that will give a suitable dose plan (i.e. one that meets 

14 all dosimetric criteria). A dose plan is computed using these parameters, and if any dosimetric criteria 

15 are not met, the process is repeated until a suitable plan is found. In this way, the dose distribution is 

16 controlled by abstract parameters. The purpose of this study is to develop an inverse planning method for 

17 HDR brachytherapy that directly controls the dose distribution through dosimetric indices. 

is We formulate HDR dose planning as an optimization model known as a mixed integer program. A mixed 

is integer program is an optimization model where some or all of the variables are restricted to taking on 

20 integer values. Dosimetric criteria can be directly incorporated into our model using integer programming 

21 constraints. Since the constraints in mixed integer programs are hard constraints, our model is guaranteed 

22 to meet the given dosimetric criteria if physically possible; it will be declared infcasiblc otherwise. However, 

23 many integer programs are difficult to solve in a reasonable amount of time [10] . This is the case for our 

24 model, and in our initial tests, reasonably sized instances of our model did not provide an optimal solution, 

25 even given several hours. 

26 HDR brachytherapy is an outpatient procedure. Catheters are inserted, a dose plan computed and 

27 delivered, and the catheters removed in the course of about one day. As a consequence, an HDR brachytherapy 
2s planning system must be able to compute a dose plan within several minutes to not interrupt clinical workflow. 

29 Since our model does not meet this runtime restriction, we develop a heuristic that does. 

30 Although our heuristic converges to a solution quickly, it has three major limitations. First, as in all 

31 heuristic approaches, our heuristic is not guaranteed to compute an optimal solution for our model. Second, 

32 our heuristic is only guaranteed to comply with certain types of dosimetric criteria. Finally, within dosimetric 

33 criteria, our heuristic tends overdose healthy tissue, even for only small improvements to tumor dosing. 

34 Despite these limitations, we show that our heuristic approach can compute dose plans quickly, and for 

35 our patient cases, did not require any iterations beyond the first. Also, contrary to current approaches, 

36 if all dosimetric criteria are not met, physicians can iterate using different dosimetric criteria rather than 
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computational parameters. A comparison of the workflow between the current approach and our approach 
is shown in Figure [T] 
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Figure 1: Clinical workflow for current approach (left) and our approach (right). With current approaches, 
physicians iterate through dosimetric criteria and parameters. In our approach, iterations are only done using 
dosimetric criteria. 

Although our model can be used for any brachytherapy cancer site where dosimetric criteria is known or 
desired, for this study we focus on the treatment of prostate cancer. In the following section, we give some 
background on inverse planning and integer programs. 



2 Background 

In HDR brachytherapy, the tumor site is referred to as the clinical target volume (CTV), and surrounding 
healthy organs are referred to as organs at risk (OAR). A physician prescribes a dose, Rx [Gy], that should 
be delivered to as much of the CTV as possible, subject to constraints on the OAR. The dose delivered to an 
entire organ is quantified by dosimetric indices. For example, the VfoQ State is the prostate volume receiving 
at least 100% Rx. 

In inverse planning, organ volume is discretized into dose points. The dose received at the dose point is 
considered representative of the dose received in the near vicinity of the dose point. If dose points are evenly 
spaced, then a dosimetric index can be computed as the sum of all the dose points receiving over the dose 
threshold times the volume that a single dose point represents. 



52 Clinical experience has determined ranges for dosimetric index values that are correlated with acceptable 

53 biological outcomes such as non-recurrence and minimal side effects to OAR. These dosimetric indices and 

54 their acceptable ranges constitute dosimetric criteria. For example, according to the RTOG-0321 dosimetric 

55 protocol [5], when treating prostate cancer, the VfQQ Statc (i.e. CTV coverage) should be greater than 90% 
se of the prostate volume and the V7 5 cctum should be less than 1 cm 3 . The objective of HDR brachytherapy 
57 planning is to determine a set of source dwell times that achieves the given dosimetric criteria. 

5s In iterative approaches to HDR brachytherapy planning such as in , suitable dose plans are computed 

59 by assigning penalties to dose points and finding a dose plan that minimizes the total penalty. For any given 

60 dose plan, dose points are penalized for receiving dose in an undesirable range. For example, a dose point in 

61 the CTV could be penalized for receiving less than Rx, and a dose point in the rectum could be penalized 

62 for receiving more than 75% Rx. 

63 The importance of individual dosimetric criterion are represented by the relative magnitude of the penalty 

64 weights. For instance, if the penalty for the CTV is three times that of the rectum, then achieving CTV 

65 coverage is more important to the physician than controlling the dose to the rectum. Penalties are useful 

66 when a standard set of dosimetric criteria is not known because it allows physicians to balance CTV coverage 

67 versus OAR exposure in a single objective. 

68 However, using penalties to achieve a prespecified set of dosimetric criteria can be difficult because penal- 

69 ties do not provide direct control over dosimetric indices. As a consequence, achieving dosimetric criteria 

70 with penalties becomes an iterative process. A set of penalties are hypothesized that will achieve dosimetric 

71 criteria; then a dose plan and dosimetry are computed. The process is repeated until a dose plan meeting 

72 all dosimetric criteria is achieved. Although a suitable dose plan can usually be found, it can take several 

73 iterations and is a time consuming process. Since standard dosimetric criteria are known for some cancer 

74 sites such as for prostate cancer, it is natural to utilize models that directly incorporate them. 

75 The model and heuristic presented in this study utilize optimization models from mathematical program- 

76 ming. In the following paragraphs we give a brief overview of these models, and describe how they relate to 

77 brachytherapy dose planning. 

78 A mixed integer program (MIP) is an optimization model of the form: 
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(MIP)Maximize c T x 
Subject to : 
Ax < b, 

Xi £ R,i £ {1, ...,k}, 
Xi £ Z, i £ {k + 1, n}, 

where the known parameters c, A, and b are nxl, mxn, and mxl matrices of real numbers, respectively, and 
x is an n x 1 vector of unknown variables. A feasible solution is an x that satisfies every constraint. The 
term c T x is called the objective function. The goal of a mixed integer program solver is to find an optimal 
solution - a feasible solution with the highest objective function value. If there are no variables in (MIP) 
that are restricted to be integers, then (MIP) reduces to a linear program (LP). LPs are efficiently solvable 
\V2\ . even on a personal computer. 

Often, the mathematical structure of MIPs makes it difficult to find a provably optimal solution for them 
within a reasonable time [lOj . As a consequence, fast heuristic, or suboptimal, approaches are required to 
make usable applications. 

Integer programs have been used to model prostate permanent-seed (PPI) brachytherapy and external 
beam dose planning [T31 HH [HI HH El [TS]. In particular, Ferris et al., 2002, used integer programs to 
constrain dosimetric indices. Our model extends these constraints to HDR brachytherapy which uses similar 
dosimetric indices. In many of these studies, extensive work was done to formulate customized algorithms to 
solve these models in a reasonable time. This study also presents a customized heuristic approach to compute 
fast but suboptimal solutions from our model. We present the model and heuristic in the following sections. 

3 Method and Materials 
3.1 Model Formulation 

We develop a mixed integer program for HDR brachytherapy planning. A description of subscripts, param- 
eters, and variables in this model are given for reference in Table [T] 

Dosimetric indices are estimated by sampling from a uniform grid of dose points that have been generated 
from evenly spaced anatomy image cross sections. The formulation of our model assumes the same set of 
dose points. The subscript "s" denotes organs, and "i" and "j" denote dose points and dwell positions, 



Term 


Description 


s 


Subscript for organs. 


i 


Subscript for dose points. 


3 


Subscript for dwell positions. 


G s 


The set of dose points in organ s. 




I he 3-D coordinates of dose point i in G 8 . 


N s 


Ine number ol dose points m G s . 




ryn ' t i ~\ i ■ , pi li ■ j • 

1 ne o-D coordinates of dwell position j . 


N T 


Ine number ol dwell positions lor tins patient. 


T) sij 


ml l i c m t i "i 

ine dose rate irom lj to P S i. 


h 


Dwell time at Tj. 




Dose at P S i. 


R s 


rril 1111 i" /^T 

ihrcshold dose for G s . 


M s 


Maximum dose for G s . 




fndicator variable for P S j. 


V s 


Dosimetric index for G s . 


L s 


Lower bound for v s . 


U s 


Upper bound for v s . 



Table 1: IPIP term definitions. 



101 respectively. The set of dose points in an organ is G s , and P S i is the 3 dimensional coordinate of a dose point 

102 in G s . The number of dose points in an organ is N s . Dwell positions are denoted by Tj , and the number of 

103 dwell positions is N T . 

104 The dose rate parameters, D si j [cGy/s], are the dose received at P si for every second the source remains 

105 at Tj. R s [cGy] is the dose required for P si to be counted in the dosimetric index for G s . The maximum 
we allowed dose for dose points in G s is M s [cGy] . L s and U s are the upper and lower bounds for the dosimetric 
10? index for G s , respectively. 

los The optimization variables are the dwell times, tj [seconds], the dose at each dose point d S i [cGy], the 

109 dosimetric indices for each organ, v s , and the indicator variables, x S i [1]. The dwell times arc continuous 

no variables that represent the source time spent at Tj. The total dose received at P S i is the sum of the 

in contributions from every dwell position, d S i = J^j^i Dsijtj. The indicator variables are binary variables that 

U2 should have the following behavior: 

J 1 if d st > IL 

x si \ 

otherwise 

N 3 

in Finally the dosimetric index for G s is the sum of all the indicator variables, v s = ^ x si . 

i=l 

U4 We dub the model Inverse Planning by Integer Program (IPIP). It is presented in its entirety below. 
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V*,ieG„(i) 

Vs,ieG s ,(2) 
V^ G G s ,(3) 

Vs,(4) 
Vj,(5) 
V*,i G G s .(6) 

115 
116 

117 The objective of IPIP is to maximize w , CTV coverage. The first constraint (1) integrates the dose at 

us P si from every dwell position. The purpose of (2) and (6) is to enforce the relationship between d si and 

us x si . For a given dose d S i at P sil if x si = 1, (2) reduces to R s < d si < M s — e, and if x si — 0, it reduces to 

120 < d si < R s — e. Conversely, if d si > i? s , then x si = 1, and if d si < R s , then x si = 0. Therefore, x si = 1 if 

121 and only if d S i > R s , and the indicator variables behave as desired. 

122 The parameter e is a small number that is included so that the value of x S i has no ambiguity when d S i = R s . 

123 More specifically, if d S i equals R s , then 1 or is valid for x S i even though it should take the value 1 according 

124 to our definition. For very small e, (2) becomes an approximation to R s x S i < d S i < R s + (M s — R s )x S i, which 

125 ensures that x S i = 1 when d S i — R s . If this level of precision is not required, then e can be omitted from the 

126 formulation. 

12? Constraint (3) integrates the dosimetric indices over all the indicator variables. Constraint (4) enforces 

128 dosimetric criteria. Since dose points lie on a uniform grid, each dose point is representative of an equal organ 

129 volume. Therefore, constraining the amount of organ volume receiving R s dose can be made equivalent to 
no constraining the number dose points receiving R s . For example if each dose point represents 0.1 cm 3 of organ 
ni volume, and no more than 1 cm 3 of organ can receive over the threshold dose, then the number of dose points 
132 in the organ that can receive over the threshold dose is 10 dose points. 



(IPIP) 
Maximize vo 
Subject to : 

N T 
3 = 1 

R s x si < d sl < R s + (M s - R s )x sl - e, 

i=l 

L„ < J, < U, 



tj>0 



x si G {0, 1} 
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133 Constraint (5) restricts the dwell times to be non-negative and (6) constrains the indicator variables to 

134 be binary. 

135 IPIP is a direct approach to finding dose plans with dosimetric criteria. However, our initial tests showed 

136 that it could not compute an optimal (or feasible solution) reliably for reasonably sized instances, even given 
is? several hours. Therefore in the following section, we develop a fast heuristic method for computing feasibly 
us solutions. 

no 3.2 IPIP Heuristic formulation 

«o The computational difficulty in solving IPIP comes from the binary variables, x si , that make IPIP a MIP. 

«i We develop a heuristic to quickly determine values for the indicator variables in IPIP. This is accomplished 

«2 by first relaxing IPIP using fewer and weaker constraints. The resulting optimization model is an LP, which 

143 can be solved quickly. The dose plan computed from this relaxation is analyzed, and based on this analysis, 

144 all x si from OAR (which have a dose upper limit) are set to 1 or 0. A doseplan is computed from the model 

145 with the set indicator variables. This dose plan is feasible for all upperbound dosimetric constraints (i.e. 
we v s < Us), but not necessarily feasible for lowerbound constraints (i.e. L s < v s ). 

14? First, the binary restriction for the indicator variables in IPIP is relaxed to allow all values between 

us and 1 (i.e. the constraint x S i € {0, 1} becomes < x si < 1). Also constraints on dosimetric index values 

149 are removed. The resulting optimization problem is the following linear program referred to as the heuristic 

150 relaxation (HR). 



N 

(HR) Maximize ^ x ai 
i=i 

Subject to : 

N T 

^^D si jtj > R s x M , Vi e G , 

N T 

tj > Vj, 
< x M < 1, Mi e G . 



151 Since dosimetric indices are not constrained, the indicator variables that make up these dosimetric indices 

152 are omitted in HR. The dose at any dose point is only limited by the maximum dose, M s . The indicator 
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153 variables for the CTV (i.e. x oi ) are retained and maximized in the objective. 

154 In general, a dose plan computed from HR will not satisfy the dosimetric criteria for the OAR in IPIP 

155 (i.e. v s > U s ). In IPIP, the value of the indicator variable x si determines if the dose at P si should be less 

156 than R s (x si = 0) or M s (x si = 1). Therefore, setting x si to or 1 is equivalent to setting the dose upper 
15? limit for P si to R s or M s respectively. In the current solution to HR, all the OAR indicator variables have 
i5s been set to 1, making the dose upper limit to these dose points M s . The next step of this heuristic is to set 
159 all but U s of the OAR indicator variables to 0. Then, v s is guaranteed to be less than U s for these organs, 
mo In this heuristic, the U s dose points in G s receiving the most dose retain a dose upper limit of M s . 

161 The remaining dose points are restricted to receive less than R s . This can be accomplished by adding the 

N T 

162 constraint ^ D si jtj < R s — e to the N s — U s dose points in G s that are receiving the least dose in the 

163 HR dose plan. The reason this allocation of indicator variables was chosen is as follows. The HR dose 

164 plan will have higher coverage than the IPIP optimal solution because it has fewer constraints. We impose 

165 additional constraints on HR so that it will be feasible for IPIP, and these additional constraints will reduce 

166 the coverage of the HR dose plan. Further restricting the coldest dose points will impact CTV coverage the 

167 least because they must be altered the least to be excluded from the dosimetric index for that organ. After 

168 these constraints are added, HR with the additional constraints is resolved as an LP. 

169 A common dosimetric criterion is for CTV coverage to be more than a certain percentage of the CTV. 

170 This is represented in IPIP by the constraint L < vq- This constraint was removed in HR and is not 

171 enforced by the additional constraints added later on. As a consequence, meeting lower bound dosimetric 

172 index constraints is not guaranteed by this heuristic. However, usually the only dosimetric index with a 

173 lower bound constraint is for CTV coverage, and this dosimetric index is maximized the objective of HR. 

174 Therefore, if the upper bound constraints are not too stringent, this criterion is likely to be met or come 

175 close. 

176 Our IPIP heuristic is summarized as follows: 

177 (1) Solve HR. 

178 (2) For each G s except the CTV, let P si * denote the N s — U s dose points receiving the least dose in G s in 

179 the HR dose plan solution. 

lso (3) For every P si e P si ^ add the constraint ^ D si jtj < R s — e to HR. 

i=i 

lsi (4) resolve updated HR to get dose plan. 

182 For the sake of brevity, we will not distinguish between IPIP and this heuristic for the remainder of the 

183 paper. 
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184 3.3 Patient Data Sets 

is5 We applied IPIP retrospectively to 20 prostate cancer patient cases. These patients were chosen to have a 

186 wide range of prostate volumes ranging from 23 to 103 cm 3 . For these patients, the physician implanted 14 

187 to 18 catheters in the prostate with transrectal ultrasound (TRUS) guidance while the patient was under 
las epidural anesthesia. Then Flexi-guide catheters (Best Industries, Inc., Flexi-needles, 283-25 (FL153-15NG)), 
189 which are 1.98-mm-diameter hollow plastic needles through which the radioactive source move, were inserted 
mo transperineally by following the tip of the catheter from the apex of the prostate to the base of the prostate 

191 using ultrasound and a stepper. A Foley catheter was inserted to help visualize the urethra. 

192 After catheter implantation, a treatment planning pelvic CT scan was obtained for each patient. Three- 

193 millimeter-thick CT slices were collected using a spiral CT. The CTV and OAR (urethra, rectum, and bladder) 

194 were contoured using the Nucletron Plato Version 14.2.6 (Nucletron B.V., Veenendaal, The Netherlands). 

195 The CTV included only the prostate and no margin was added. When segmenting the bladder and rectum, 

196 the outermost mucosa surface was contoured. The urethra was defined by the outer surface of the Foley 

197 catheter, and only the urethral volume within the CTV was contoured. The OAR were contoured on all CT 

198 slices containing the CTV and at least two additional slices above and below. Implanted catheters were also 

199 digitized. 

200 For the contoured anatomical structures, dose points were generated by sampling from a uniform grid 

201 with 2 mm spacing in the x-y direction and 3 mm spacing in the z direction. Dose points were also generated 

202 in the body tissue (space between organs) with 4 mm spacing in the x-y direction. For the 20 cases, the 

203 total number of dose points for IPIP ranged from 8860 to 25288. To reduce the computation required, every 

204 dose point more than 3.5 cm from the xy-centroid of the dwell positions were omitted in the optimization. 

205 However, they were included when computing the dosimetric indices. Our tests showed that removing these 

206 dose points had no effect on the dosimetric index values. The number of dose points after removing dose 

207 points ranged from 6144 to 15568. 

208 3.4 Dose rate calculations and clinical criteria 

209 The dose-rate contribution to a dose point from a source dwell position is a function of the distance between 

210 them. The dose-rate parameters were calculated as specified in the AAPM TG-43 dosimetry protocol [2H] , 
2u The radioactive material used in the source was 192 Ir, and the prescription dose was 9.5 Gy. 

212 The clinical criteria used in this study can be found in Table [5J The specifications for the V^oo statc , 

213 Vgf hra , V^ 5 octum , and V^ 5 laddcr are defined by RTOG-0321 pj. RTOG specifies that the Vgf 411 ™ be much 

214 less than 1cm 3 . We interpreted this to mean less than 0.1cm 3 . 
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Ind6x 


R,g qui rc nic nt 


V 100 


> 90% 


"v tJt'i'Os t ;| | p 

V 150 


< 45% 


■\ j U rctlira 
V 120 


< 0.1 cm 3 


\j Urethra 
V 150 


= cm 


■frHectum 
v 75 


< 1 cm 3 


"\rHcctum 
v 100 


cm 3 


\tH ladder 
V 75 


< 1 cm 3 


\tH ladder 
v 100 


= cm 3 


v 200 


= cm 3 



Table 2: Dosimetric Criteria. 

21 5 The Vf 5 r £ stato is not explicitly constrained by the RTOG-0321 protocol. The Vf 5 r £ statc is restricted by the 

216 homogeneity index (HI) where: 



HI 



V Prostate \rProstatc 
100 v 150 



"\ /Prostate 
v 100 



2is It is generally preferred that HI > 0.6; however, lower values of HI are acceptable if they allow for higher 

219 CTV coverage. For this study, we constrain HI to be greater than 50% to maintain some control over the 

220 Vf 5 r Q Statc while not being overly restrictive. Since we expect target coverage over 90%, this restriction on HI 

221 can be enforced with V^ stat0 < 45%. 

The restrictions that the V^ thra , V^ ctum , Vf(] addor , and Vf ° Q dy be equal to are not specified by RTOG 

223 but are considered preferable when possible at our clinic. The preference of the V^f 5 ' comes from the desire 

224 to keep hot spots localized within the CTV. 

225 These dosimetric index constraints represent our dosimetric criteria. They are summarized in Table [2] 

226 The parameters used for IPIP to reflect these dosimetric criteria can be found in Table [3| The values 

227 of 8 and 83 for U s represent the number of dose points in 0.1 cm 3 and 1.0 cm 3 , respectively, based on our 
22s grid spacing. The value of U s for the Vf5Q StatG is 45% the number of dose points in the prostate. The dose 

229 to prostate dose points should be unrestricted so we have used an unrestrictively high number, 20000 cGy. 

230 We did this to avoid using infinity, which creates numerical problems with our optimization solver. The 

231 dosimetry in our results shows no dose points receiving this dose level for any patient. 



232 3.5 Method evaluation 

233 We used Matlab v. R2008b (Mathworks Inc.) to compute dose plans from IPIP. The linear programming 

234 optimization was done using the Matlab interface for the Mosek Optimization Toolbox v. 5, a medium to 

235 large scale optimization package. All computations were performed on a personal computer with Intel(R) 
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Prostate 
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950 


2000U 


iVo 


Prostate 


1 


1425 


20000 


.45 N 


Urethra 


2 


1140 


1425 


8 


Rectum 


3 


712 


950 


83 


Bladder 


4 


712 


950 


83 


Body 


5 


1900 


1900 






Table 3: IPIP Parameters 





Approach 


Dosimetric Index 


IPSA 


IPIP 


wProstatc 
v 100 


100% 


100% 


~\ /Prostate 
v 150 


100% 


100% 


t rUrcthra 
v 125 


85% 


100% 


t r Urethra 
V 150 


50% 


100% 


^ rlicctum 
V 75 


80% 


100% 


"\ rlicctum 
v 100 


55% 


100% 


"\rBladdcr 
V 75 


65% 


100% 


"\rBladdcr 
v 100 


100% 


100% 


v Bod y 

v 200 


5% 


100% 


HI 


100% 


100% 


All Requirements 


0% 


100% 



Table 4: IPSA and IPIP compliance with dosimetric criteria. Percentage of patient dose plans (out of 20) 
that complied with each constraint for each approach. Note that the dosimetric constraints supplied were 
hard constraints, meaning even the slightest deviation from the limit is considered a failure. 

236 Core™ 2 Duo CPU 1.67 GHz processor, 3 Gib RAM, and the Windows 32-bit operating system. We record 

237 the compliance of IP2H dose plans with our dosimetric criteria and the running time. 



4 Results 

239 IPIP satisfied all our dosimetric requirements on the first iteration. For comparison, we also computed 

240 dose plans from Inverse Planning Simulated Annealing (IPSA), a clinically deployed dose planning system 

241 used at our clinic. All computational parameters for IPSA were the same as in Alterovitz, 2006 [19 . No 

242 computational parameters from IPSA were manipulated after the first iteration. 

243 The compliance rate (out of 20 patients) for each individual dose objective is summarized in Table |4j 

244 Note that our constraints on dosimetric indices are hard constraints so even the slightest deviation from the 

245 requirement is considered a failure. The average CTV coverage from IPIP was 95%. 

246 On average, IPIP computed a dose plan within 30.1 seconds. The maximum runtime was 86 seconds. Our 

247 IPIP heuristic is made up of (1) solving a linear program, (2) sorting doses within each OAR, and (3) solving 

248 another linear program. Each one of these steps is polynomial time solvable, meaning that the number of 
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X^/UolJ.lJ.Lv \j±. IVj 111U.I-/A 




"\ j Frost ate 
v 100 


78% 


\ j Frost ate 
V 150 


33% 


t r Urethra 
V 125 


0.01cm 3 


"\ /" Hoc turn 
v 75 


cm 3 


t Bladder 
V 75 


cm 3 



Table 5: Set of dosimetric criteria that cannot be met with IPIP. The y^g * 1 ™ and v!? 5 were required to 
be 0. For this case, CTV coverage was only 78% which does not meet our 90% requirement. Less stringent 
dosimetric criteria must be given in the next iteration. 



Dosimetry 


Dosimetric Criteria 1 


Dosimetric Criteria 2 


"\ rProstatc 
v 100 


96% 


91% 


"\/Prostatc 
v 150 


33% 


30% 


"\ 7" U ret lira, 
v 125 


0.01cm 3 


cm 3 


t rllcctum 
V 75 


0.92 cm 3 


0.45 cm 3 


"\rBladdcr 
V 75 


0.95 cm 3 


0.44 cm 3 



Table 6: Comparison of different dosimetric criteria for a single patient. For Dosimetric Criteria 1, dosimetric 
criteria was the same as in our main experiment. For Dosimetric Criteria 2, the y^ ctum and v B 5 laddcr are 
restricted to be less than 0.5 cm 3 , which is more stringent. By restricting the V^ 5 cctmn and v B 5 laddor , CTV 
coverage drops by 5% but with large reduction to OAR exposure. 

249 computations grows slowly (i.e. as a polynomial) with respect to the size of the input. As a consequence, 

250 the runtime performance of our IPIP heuristic is reliable for a more general audience of patients given that 

251 they have a similar number of dose points. 

252 For comparison, IPSA computed a dose plan within 5 seconds on average, with a maximum runtime 

253 of 9 seconds. However, IPSA would have required further iterations to achieve dosimetric criteria. These 

254 additional iterations are not included in the running time results for IPSA. 

255 We also conducted side tests to demonstrate cases in which IPIP docs not behave as desired. First 

256 we showed a set of dosimetric criteria that cannot be met with IPIP by insisting that the y^,? * 1 ™, and 

257 v^ 5 cctum be cm 3 . Otherwise, dosimetric criteria were the same. The results for a single patient are shown 

258 m Table[5| CTV coverage, which is lower bound constrained, is only 78% which does not meet the 90% 

259 requirement. However, all other criteria are fulfilled. In this case, additional iterations would be required 

260 using less stringent dosimetric criteria. 

261 Second, we showed a case where a dose plan computed from IPIP complies with dosimetric criteria but 

262 may not be what the physician desires. We computed two dose plans from IPIP using different sets of 

263 dosimetric criteria. The first dosimetric criteria set was the same that was used in our main experiment. 

264 The second dosimetric criteria set was more stringent with v^ 5 cctum and v B 5 laddcr less than 0.5 cm 3 . All other 

265 dosimetric criteria were the same as for the main experiment. The results are summarized in Table [6] 

266 CTV coverage for the first dosimetric criteria is 96%. To achieve this coverage, almost all of the dose 
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267 allowance to the rectum and bladder are utilized (i.e. v^ 5 cctum and V75 addor is close to 1 cm 3 ). With the 

268 second dosimetric criteria, CTV coverage drops to 91%, but the overdosed volume to the rectum and bladder 

269 is cut in half. 

270 Since the objective of IPIP is to maximize coverage within limits to OAR dose, IPIP will tend to utilize 

271 all of the OAR dose available to it. More specifically, IPIP will make large increases to OAR dose, within 

272 its limits, even to make only small improvements to CTV coverage. In this case, the drop in CTV coverage 

273 (5%) between the two dose plans is worth considering, but a physician may still prefer the second dose plan. 

274 However, determining the second dose plan cannot be accomplished without iterating through alternative 

275 dosimetric criteria. 

276 5 Conclusion 

277 IPIP is a new approach that is more intuitive for physicians and less reliant on manual fine tuning. IPIP 

278 allows physicians to directly specify desired dosimetric indices, rather than system-specific computational 

279 parameters. Our results demonstrate that IPIP quickly generates dosimetry plans that are consistent with 

280 RTOG-0321 standard dosimetric criteria. Further study is required to determine if IPIP can produce similar 

281 performance for a more general group of patients and dosimetric criteria, including other cancer sites such 

282 as GYN. 
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